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ABSTRACT 

The nearest-neighbour XY spin glass on a hypercubic lattice in four 
dimensions is studied by Monte Carlo simulations. A finite-size scaling 
analysis of the data leads to a finite temperature spin glass transition at 
T c = 0.95 ± 0.15. The critical exponents are estimated to be v sg = 0.70 ± 
0.10 and rj sg = —0.28 ± 0.38. The results imply that the lower critical 
dimensionality for the XY spin glass is less than four. 



2 



Recently there has been considerable interest in the behaviour of short- 
range vector spin glasses [1-8]. Although there is now convincing evidence 
[2-4] that XY spin glasses exhibit only a zero-temperature phase transition 
for d = 2 and 3, the location of the lower critical dimension, g^, remains 
controversial. A zero-temperature study by Morris et al [2] suggests that 
dL = 4. The validity of Nishimori and Ozeki's [5] attempt at a Mermin [9] 
type argument has been questioned [6] and it is claimed by Schwartz and 
Young [7] that all one can actually show is that d^>2. 

A recent Migdal-Kadanoff renormalisation-group analysis [8], on the 
other hand, claims that the XY spin glass orders at a finite temperature for 
d = 4 and, therefore, < 4. The XY spin glass, unlike the corresponding 
gauge glass, possesses a 'reflection' symmetry. It has therefore been argued 
[10, 11] that the two models belong to different universality classes. Recent 
computer simulations [11] of the gauge glass in four dimensions clearly show 
a finite temperature transition. 

As far as we are aware, to-date the XY spin glass has not been studied 
by simulations for d = 4. We attempt to fill this gap by presenting in this 
letter the first results of extensive Monte Carlo simulations of the four di- 
mensional XY spin glass with ± J nearest neighbour interactions. Applying 
finite-size scaling theory [12], we shall find evidence for a finite-temperature 
glass transition. 
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The Hamiltonian for the model simulated is given by 

<ij> <i,j> 

with < 0j < 27T for all i. The planar spins, Sj(= are situated 

on every site of a four dimensional hypercubic lattice of size L 4 (L = 2,4 
and 6). The summation runs over nearest-neighbour pairs only and the 
interactions, Jij, are independent random variables selected from a binary 
±1 distribution. As usual, the temperature is given in units of the nearest 
neighbour interaction. We work with full periodic boundary conditions. 

In the simulations we study the spin-glass susceptibility, Xsg, which is 
defined by 

= Hf, (2) 

where < ... >t indicates a thermal average, [...]j denotes an average over 
the disorder, N = L 4 and 

4f = £[< 1% >t]j- (3) 

In equation (3) the tensor variable, q^, is defined in terms of the overlap 
between two replicas 1 and 2, 

= ~^Yl S l» S lv (v,v = x,y). (4) 

i 
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Higher order correlations such as q^Q can also be written in terms of q^, 



namely 

,(4) = 



Vs7 = I< qIv&p >t]j- (5) 



Rather than use q^g and q^ g \ it is far more convenient to work with the 
dimensionless Binder parameter [12] defined by 

a (4) 

According to finite-size scaling theory [12], near T c we expect the Binder 
parameter to scale as 

g sg (L,T)=g sg (L 1 ^(T-T c )), (7) 

where v sg is the correlation length exponent and g sg is a scaling function 
satisfying 



g sg (L 1 ^(T-T c )) = {° 1 



forT>T c ,L^oo () 

forT<T c ,L^oo ^ 



provided that the ground state is non-degenerate. Whereas for a finite- 
temperature spin glass transition plots of g sg (L,T) versus T for different L 
should intersect at T c , for a transition at zero temperature we expect the 
curves to meet each other at T = only. 

The scaling form for the spin glass susceptibility is given by 

Xsg {L,T) =L 2 -^x sg (L 1 ^(T-T c )), (9) 



where the exponent t] sg describes the power-law decay of correlations at the 
transition temperature and x sg is now a scaling function. It follows from 
equation (9) that 

Xsg (L,T c )~L 2 -^. (10) 

As the co-ordination number for our model is Z = 8, the mean- field values 
of T c and the exponents mentioned above are 

T c m/ « = y/2, v mf = 1/2 and r] mf = 0. (11) 

We now describe our Monte Carlo simulations and discuss the results. Dur- 
ing the simulations, which were performed using the conventional Metropo- 
lis [13] technique on a Cray Y-MP8 and a J90, we actually work with 
discrete spins. For technical reasons, the spins were discretized to occupy 
256 equally spaced orientations in the plane. Furthermore, a variant of 
multispin coding [14] was used to store 7 (discrete) spins in one word and 
the lattice was composed of two inter-penetrating sub-lattices. As a conse- 
quence, each update of the lattice allows us to update 14 separate samples 
(or, alternatively, 7 pairs of samples) at the same time. 

We follow Bhatt and Young [12] and compare the spin glass correlations 
obtained from 2 independent replicas with the same set of bonds with those 
from a single replica at 2 different times. Equilibrium is assumed only if 
the values agree within the statistical error. 

The number of Monte Carlo steps, tq(L,T), required to achieve equi- 
librium depends on both the system size L and the temperature T. To(L, T) 
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sets upper and lower limits on the values of L and T, respectively, that can 
be studied. In our simulations we found that To (2, 0.4) « 1000, To (4, 0.7) « 
5000 and To (6, 0.75) « 13000 sweeps; equilibration problems prevented us 
from going to lower temeperatures. The number of bond configurations we 
generated to average over the disorder also varied with L. Typically, we 
considered 7000 (L = 2), 250 - 2300 (L = 4), 100 - 500 (L = 6) pairs of sam- 
ples for each temperature considered. (However, note that for the lowest 
2 temperatures (T = 0.775 and 0.75) for L = 6 we averaged over only 56 
pairs of samples in each case.) In total, the simulations presented in this 
work took approximately 600 hours of CPU time on the two supercomputers 
mentioned above. 

In figure 1 we plot Xsg against T for the 3 different values of L consid- 
ered. The statistical error-bars were evaluated from the sample-to-sample 
fluctuations and are only displayed in those cases where they exceed the 
size of the points. 

The Binder parameter is plotted against the temperature in figure 2. 
Although the curves clearly intersect at a finite temperature, the point of 
intersection is not unique. This is probably due to corrections to finite-size 
scaling and the statistical error in g sg . As we have data for 3 values of L, 
we obtain 3 intersection temperatures, T c Ll ' L2 , where L\,L<i = 2,4,6 and 
L\ 7^ L/2- For the data presented in figure 2 we note that there is a small 
downward shift in the value of T^ 1,Ll2 for increasing L\ and L^. Clearly, 
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to establish whether the shift is significant or not, it is highly desirable to 
obtain additional data for the Binder parameter for larger lattices and lower 
temperatures. However, we note that there is some evidence of a finite- 
temperature transition as the curves clearly splay out below the intersection 
point. From figure 2 we estimate the spin-glass transition temperature to 
be T c = 0.95 ± 0.15. Our value of T c (« 0.7 T c m/ ) agrees well with the value 
of T c « 0.9 obtained recently by Nobre et al [8]. It is also surprisingly close 
to the transition temperature found by Reger and Young [11] for the four 
dimensional gauge glass. 

Having obtained an estimate for T c , a log-log plot of Xsg(L, T c ) against 
L is expected from equation (10) to have a slope of 2 — rj sg . Our results are 
consistent with this but the uncertainty in T c leads to a large error in rj sg 
and we estimate that 2 — rj sg = 2.28 ± 0.38. 

To fix the second independent exponent, v sgi we display in figure 3 
a scaling plot of x sg (L,T)/L 2 - 28 against {L X / V ^{T - 0.95)). To see how 
sensitive the scaling plot is to the value of v sg , we have tried a range of 
values. As a result we estimate v sg = 0.70 ± 0.10. As can be seen from 
figure 3, the data scale reasonably well for v sg = 0.70. Once again, our 
value for the correlation exponent is remarkably close to the one found for 
the gauge glass in 4d [11]. 

Finally, in figure 4 we show a scaling plot of the data for the Binder 
parameter for the same values of T c and v sg as the ones used in the plot for 
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figure 3. We see that the data for g sg do not scale as well as those for Xsg- 
The quality of the data collapse does not improve for other possible values 
of T c and v sg . (A correction to scaling as suggested by Bokil and Young [15] 
also fails to make any significant difference to the plot.) A similar behaviour 
was found by Kawamura [4] in the three-dimensional XY spin glass. 

It has been assumed that the XY spin glass and the gauge glass belong 
to different universality classes as the latter does not share the reflection 
symmetry of the former. As noted, our results are remarkable in their 
similarity to those found earlier by Reger and Young [11] for the vortex 
glass in four dimensions. It is felt that this unexpected feature requires 
further investigation. 

To conclude, we have presented numerical evidence that the XY spin 
glass has a finite temperature glass transition in four dimensions. We have 
estimated the critical temperature and the critical exponents. Further work 
is required to confirm the transition temperature and obtain more accurate 
values for the exponents. Our results are in agreement with the analytic 
approximation carried out by Nobre et al [8]. They are also surprisingly 
similar to those found earlier by Reger and Young [11] for the four dimen- 
sional gauge glass and imply that cIl < 4 for the XY spin glass. 

Work is underway to investigate the chiral-glass [4, 15] behaviour of 
the model using the vortex representation [15]. 

The simulations were performed on a Cray YMP and a J90 at the 
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Rutherford Appleton Laboratory through an Engineering and Physical Sci- 
ences Research Council (EPSRC) research grant (Ref: GR/K/00813). 
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FIGURE CAPTIONS 



Figure 1 

A plot of the spin glass susceptibility, Xsg, against the temperature for 
L = 2, 4 and 6. The lines are just guides to the eye. 

Figure 2 

A plot of the Binder parameter defined in equation (6) against the temper- 
ature for L = 2, 4 and 6. The lines are just guides to the eye. 

Figure 3 

A scaling plot of Xsg/L 2 ' 7139 versus L^^^iT - T c ) with r) sg = -0.28, v sg = 
0.70 and T c = 0.95. See equation (9) in the text. 

Figure 4 

A scaling plot of the Binder parameter g sg versus L 1 ^ 89 ^ — T c ) for v sg = 
0.70 and T c = 0.95. 
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